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Abstract 



It is shown that in a model in which up-type and down-type fermions 
acquire mass from different Higgs doublets, the anthropic tuning of the 
■ Higgs mass parameters can explain the fact that the observed masses 

CT) ' of the d and u quarks are nearly the same with d slightly heavier. If 

Yukawa couplings are assumed not to "scan" (vary among domains), 
this would also help explain why t is much heavier than b. It is also 
pointed out that the existence of dark matter invalidates some ear- 
lier anthropic arguments against the viability of domains where the 
Standard Model Higgs has positive /it 2 , but makes other even stronger 
Q_i' arguments possible. 



1 Introduction 



The mass parameter of the Higgs field in the Standard Model (/x 2 ) gives the 
appearance of being "anthropically tuned" [TJ. That is, if one imagines the 
other parameters of the Standard Model to be fixed, and considers what 
the universe would look like for different values of /i 2 , one finds that organic 
life may only be possible if /i 2 is negative and has a magnitude very close 
to the value actually observed. From the physics point of view this might 
be just a coincidence, though a remarkable one. On the other hand, as 
pointed out by Ref. pQ, it might have a physical explanation in the context 
of "many-domain", "multiverse" , "landscape" scenarios. (These are equiv- 
alent names for the same idea. For recent reviews see [2].) An explanation 
of that sort would account for the closeness of the strong-interaction and 
weak-interaction scales (M s t rong /Mpi ~ 10 -19 and M wea k/Mpi ~ 1CP 17 ), 
something not yet explained by any other scenario. (The nearness of the 
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weak-scale to the strong scale results in the lightest quarks having masses 
small compared to Aqcd, which in turn leads to the existence of nonpertur- 
bative bound states of quarks, pseudogoldstone pions, and the the richness 
of hadronic and nuclear physics, and consequently of chemistry.) 

In this paper we extend the analysis of [1] to a slightly more general 
Higgs structure and show that a possible anthropic explanation of the fact 
that m u /md ~ 1 emerges. We also reconsider the crucial case of positive 
fj?, and point out that one of the anthropic arguments proposed in pQ to 
exclude most of the > region is invalidated by the existence of dark 
energy (which had not yet been discovered when pQ was written). However, 
we show that the existence of dark energy makes possible a very different 
set of arguments that reach even stronger conclusions. 

The papers of Ref. 1 assumed that the only parameter of the Stan- 
dard Model that varies among domains and is anthropically tuned is /i 2 . 
Several facts suggest the possibility that the observed value of /i 2 may be 
acounted for in a different way than the values of the other parameters of 
the Standard Model. First, /i 2 is the only dimensionful parameter of the 
Standard Model Lagrangian. Second, it is the most highly tuned of the 
Standard Model parameters, being 10 -34 of its "natural" value. (The next 
most tuned parameter is 9, which is less than 10~ 9 of its natural value.) 
Third, the smallness of /i 2 is so far the most intractable of the naturalness 
problems of the Standard Model. Various plausible mechanisms of a more 
conventional sort have been proposed for explaining the smallness of other 
Standard Model parameters. (For example, the smallness of 8 can be ex- 
plained by the Peccei-Quinn mechanism [3] or by spontaneously broken CP 
or P [4]; and various symmetry schemes have been proposed to explain the 
smallness of the Yukawa couplings of the light quarks and leptons). By con- 
trast, attempts to explain the smallness of the weak scale by technicolor or 
similar ideas are plagued by a variety of well-known difficulties. And low 
energy supersymmetry does not by itself explain the magnitude of the weak 
scale, though it protects it from radiative correction. 

Since the focus of Ref. 1 was on the mass parameter of the Higgs field, 
the question naturally arises how things would be different if there were 
two or more Higgs doublets. Two Higgs doublets appear in a wide variety 
of theoretical contexts: in theories with supersymmetry; in Peccei-Quinn 
models; and in many grand unified models (such as SO(W), if the 10 of 
Higgs fields is assumed to be complex). The simplest possibility, which 
is typical of the above-mentioned scenarios, and which avoids problems of 
Higgs-mediated flavor-changing processes, is that one doublet (H u ) couples 
to the up-type quarks, while the other doublet (H^) couples to the down- 
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type quarks and charged leptons [5]. It is this situation we will analyze in 
this paper. Thus, at high scales we have 

C-Yuk = (Y u uu + y c cc + Y t tt)H u 

+ (Y d dd + Y s ss + Y b bb)H d + (Y e ee + YtfH + Y T TT)H d . 

With two Higgs doublets, there are several dimensionful Higgs-mass param- 
eters that may vary among domains, with the consequence that (H u ) = v u 
and {H d ) = v d may vary independently. This means that there is the possi- 
bility of explaining not only the magnitude of the weak scale, but also the 
relative magnitudes of the masses of the up-type and down-type quarks. In 
particular, we shall see that an explanation of the size of m u /md and thus 
a partial explanation of mt/mb emerges. Much of our analysis would apply 
equally to models with or without low-energy super symmetry, and would ap- 
ply in fact to any model where in effect both v u and v d vary among domains. 
Thus, the "anthropic" considerations given in this paper are not proposed 
as an alternative to low-energy supersymmetry. Indeed, they may shed light 
eventually on supersymmetry-breaking and the fi problem. However, in this 
paper we will consider for purposes of concreteness a non-super symmetric 
two-Higgs-doublet model. 

In such a model, the parameter /x 2 is replaced by a two- by- two Hermitian 
mass matrix Mfj, i,j = u,d, with four real parameters. To make the weak- 
scale small, still only one tuning is required, namely of the determinant of 
M 4 2 , since only one Higgs doublet needs to remain light to break SU(2)l x 
U(1)y- The other doublet has no reason to be light, and in fact would 
be expected to be of superlarge mass. The doublet that is light (whose 
mass-squared we shall call /U 2 ) is a linear combination of H u and H d : H < = 
sin PH U + cos f3H d , where H d = ^H^. The mixing angle f3 enters into 
the relative magnitudes of the up-quark masses and down-quark masses. 
In particular, the ratios m u /m d and mt/mb are proportional to tan/?. If 
M 4 2 varies in the landscape (or "scans"), then so does tan (3. We shall see 
that "anthropic" considerations strongly favor the observation of values of 
m u /m d that are of order (but somewhat smaller than) 1 and thus of values 
of mt/mi, much larger than 1 (assuming, as we shall, that quark and lepton 
Yukawa couplings are fixed, i.e. do not vary among domains). 

It is noteworthy that in published models or scenarios that do not in- 
voke "anthropic" or "landscape" considerations the relations m u /m d ~ 1 
and mt/mb S> 1 are generally not accounted for but are simply assumed to 
hold. (However, there have been a few attempts to explain m u ~ m d in more 
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conventional ways, such as [6].) It seems that anthropic and non-anthropic 
explanations are somewhat complementary: those relationships that seem 
most resistant to conventional dynamical or symmetry explanations (which 
would include not only m^/m^ ~ 1, but also M wea kjMpi ~ M s t rong jMpi 
and A/M pe ~ 1(T 120 ) are also those that are most straightforwardly ac- 
counted for "anthropically" . 

In this paper, as in Ref. 1, we assume that the dimensionless parameters 
of the Lagrangian are fixed, i.e. do not vary among domains. (More pre- 
cisely, we assume that their values at some high scale are fixed, since their 
low-scale values will be slightly influenced through renormalization group 
running by the values of v u and Vd-) If this is assumed, then the inter- family 
ratios of Yukawa couplings are known, e.g. the ratios of Y u /Y t , Y c /Y t , Y^/Yb, 
and Yj/Yj,. We will also assume that the Yukawa couplings of the third fam- 
ily are all of order 1. This is plausible, given that schemes that unify all the 
fermions of a family, such as SO(W), relate these three couplings; and the 
simplest such schemes give Y t = Yj, = Y T . 

The relative values of the Yukawa couplings of the up-type and down- 
type quarks are shown in a log plot in Fig 1. 

log(Y) 
- 



Y b ,Y t 



-1 



-2 



-3 



-4 



Figure 1: The Yukawa couplings of the quarks if Yj, ~ Yt ~ 1. The up- type 
quarks have a much steeper family hierarchy. The scale on the left is log 10 Yj. 
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Since the mass hierarchy among the up-type quarks is much steeper 
than that among the down-type quarks, normalizing Yj ~ implies that 
Y d > Y u . That is, 



~ 80. 



(2) 



In the scenario we are discussing, the relative normalization of the up-type 
masses and down-type masses is free to "scan", i.e. vary among domains. 
The consequence of this, as we shall see, is that m d jm u is anthropically 
selected to be close to but somewhat larger than 1, and therefore mt/rrib 
ends up (in anthropically viable domains) being much larger than 1, as 
shown in Fig 2. It is in this sense that one can anthropically "explain" the 

log (m) 



m t 



m b 



-m d 



Figure 2: If v u /vd is anthropically set so ~ 1, then m t 3> due 

to the steeper hierarchy of the up-type quarks. The scale on the left is 
log 10 (mi/MeV). 



largeness of mt/mf,. Of course, this assumes the pattern of interfamily mass 
ratios shown in Fig 1, i.e. the greater steepness of the hierarchy for up- 
type quarks than for down-type quarks. There exist many interesting and 
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plausible theoretical schemes for explaining that pattern. One example is 
the model of [7] , where the hierarchy among the up- type quarks is quadratic 
and the hierarchy among the down-type quarks is linear in the same small 
parameters. 

In this paper, we shall call domains of the universe where organic life 
based on chemistry is possible "viable", and we shall call the values of pa- 
rameters in such domains "anthropically allowed" . In section 2 of this paper, 
we shall examine the consequences for "viability" of different values of m u 
and ma and graphically display the anthropically allowed and forbidden re- 
gions of the In m u -ln plane. In section 3, we shall look at the probability 
distributions of (m u , m^) that arise from certain simple assumptions about 
the Higgs mass matrix Mfj. In section 4 we will use the probability distri- 
butions derived in section 3 to analyze the fi 2 > region of parameter space 
in light of the existence of dark energy and the anthropic requirement that 
it be fine-tuned to be small. 

2 Physics in the lnm u — lnm^ Plane 

In this section we will discuss anthropic constraints on m u and m^. Most of 
these constraints have been discussed in slightly different contexts in pQ and 
[8], and much of the discussion in this section makes use of their results. 

The vacuum expectation v that breaks the weak interactions and gives 
u and d their masses is that of the lighter Higgs doublet H < , whose mass- 
squared parameter we are calling /i 2 . We suppose that has a su- 
perlarge "natural" scale M#, which is of order Mpi (or Mqut, it does not 
matter for our analysis which). There are two qualitatively very different 
cases, fi 2 < and /j, 2 > 0. If /j, 2 < 0, as in our domain, then v ~ a 1 / / |/j^T 
and is naturally of order M*. In this case, smaller v corresponds to more 
fine-tuning of M 4 2 . On the other hand, if /i 2 > 0, the breaking of the weak 
interactions is only induced by the coupling of H < to the quarks, which 
acquire condensates through the dynamical breaking of chiral symmetry by 
QCD, so that v ~ f^/fJ- 2 - In this case, the natural scale of v is f^/M 2 , and 
larger v corresponds to more fine-tuning. 

One sees that the maximum and minimum values of v are different by 
a factor of (M*//„-) 3 . In other words, if ~ Mp£, v can range over 60 
orders of magnitude. In fact, v u and Vd each span the same range. This 
enormous space of possible values of m u and is shown in Fig. 3. In 
that figure are also plotted the contours of equal probability density in the 
landscape (i.e. equal amount of fine-tuning of parameters), which will be 
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x u = ln(m u /M P1 ) 

Figure 3: The possible range of mj, i = u,d, is from Yif%/Mp £ to YjMp^. 
The blue dot shows the observed values. The contours are of probability 
density in the landscape, red (green) green for the [i 2 < (> 0) region. 
Neighboring contours differ by e 20 in probability density, increasing in the 
direction of the arrows. 



derived in section 3 under the assumption that the parameters of M 4 2 have 
a flat probability distribution in the landscape in the interval [— M 2 , +M 2 ]. 
The red curves are contours of equal probability density in the fi 2 < region, 
and the green curves are contours of equal probability density in the fi 2 > 
region. Neighboring contours differ in probability density by a factor of e 20 , 
with the direction of increasing probability density indicated by the arrows. 

The position of the "anthropically allowed" region in Fig. 3 in which our 
domain is located is indicated by the blue dot. In fact, the viable region is 
so small that it is smaller than the blue dot in Fig. 3. Fig. 4 is a blow-up 
of that region of parameter space. 

Since we are assuming that the Yukawa couplings of the quarks and 
leptons do not vary among domains, and that both the down-type quarks 
and charged leptons get mass from the same Higgs doublet Hj, the electron 
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mass is given by m e = mid, where r is a fixed ratio (up to logarithmic 
dependence on v u and Vd through the running of parameters). We take 
m-uO = 4 MeV and nido = 7 MeV, where the subscript zero refers to the 
observed values of these parameters, i.e. their values in our domain of the 
universe. Thus r = 0.07. 

2.1 The case of < 

First, we will discuss the case of negative fj?. If u? is negative and smaller 
in absolute value than in our domain, that involves more fine-tuning and 
is thus a less probable region of parameter space. Our interest is therefore 
primarily in fi 2 negative and greater than (in absolute value) or comparable 
to its value in our domain. That corresponds to values of m u and nid greater 
than or comparable to their observed values. 

If nid is sufficiently large, then all d quarks will beta decay into u quarks 
even inside hadrons, so that the only stable hadron will be A ++ (uuu). 
"Sufficiently large" in this case means large enough to overcome the extra 
energy required to put three quarks of the same flavor in a totally antisym- 
metric state, which energy we call -Eatv- For v u and Vd near the observed 
weak scale vq, Ean ~ 300 MeV. (However, there is a weak dependence 
of i?A7V on v u and Vd through the renormalization-group running of a s , 
due to quark thresholds.) The condition for all d quarks to beta decay is 
m<d > n^u + m e + E&n, or (1 — r)vrtd > m u + -Eaat- This is the region above 
the upper blue curve in Fig. 4, which is a log-log plot. For m u and nid 
masses very large compared to the QCD scale, this curve asymptotically ap- 
proaches the straight line (l-r)nid = m u , which is a line of slope 1 in Fig. 4. 
For m u and nid small compared to the QCD scale, the curve asymptotically 
approaches the horizontal line nid = ^Aiv/(1 — r )- 

In the region where all the d quarks decay, not only will the A ++ be the 
only stable hadron, it will surely also be the only stable nucleus. (In this 
regime, nid is large compared to the QCD scale, so that the pions are not 
pseudo-goldstone particles and the range of the nuclear potential is shorter 
than in nuclei of our own domain. Moreover, there is a strong Coulomb 
repulsion between two A ++ .) These A ++ will clothe themselves with a pair 
of electrons, and act chemically like helium. With the only element being 
helium-like, it is extremely implausible that life based on chemistry could 
exist. 

On the other hand, if m u is sufficiently large compared to m^, then all 
the u quarks will decay to d, and one will have a domain in which the only 
stable hadron and nucleus will be a A~ (ddd). These will clothe themselves 
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-54 -52 -50 -48 -46 -44 -42 -40 
x u = ln(m u /Mpi) 

Figure 4: A blowup of the region near our domain in Fig. 3. The curves are 
anthropic bounds on m u and discussed in the text. The green shaded 
region is potentially viable. The dot is our domain. 



with positrons (which will presumably be more numerous than electrons 
in such a domain, due to the beta decay of all the u quarks), and thus 
act chemically like hydrogen. An all-A~ universe is also extremely unlikely 
to give rise to chemistry-based life. The all- A - universe corresponds to 
m u > m d(^ + r ) + -Eaat? which is the region below the lower blue curve 
in Fig 4. (For large quark masses this curve it is asymptotic to the line 
(1 + r)m^ = m u , which is a line of slope 1 in Fig. 4.) 

Thus we can rule out almost all of the three quadrants in Fig. 4 where 
either m u or is large compared to its observed value. There is a thin 
ribbon between the two blue curves where both m u and irid are large com- 
pared to the QCD scale and (1 — r)m^ < m u < (1 + r)m^ where neither 
d — > u + e~ + V e nor u — > d + e + + u e can happen. However, this region is 
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unlikely to be viable, since A = 2 nuclei are not bound there, as shall be 
seen. 

More of the parameter space is ruled out as viable if we take into account 
the conditions that both neutrons and protons should be stable inside nuclei. 
If (A, Z) denotes a nucleus of mass A and atomic number Z, then the process 
(A, Z) —* (A—l,Z)+p + e~+V e is energetically allowed if — m u — m e — 
Am em > B, where B is the binding energy of the proton in the nucleus, 
and Am em is the electromagnetic contribution to the proton-neutron mass 
difference. Following PQ, we take Am em = 1.7 MeV, and we take 10 MeV as 
a typical value of B. (Of course, for very large values of v u and v^,, the values 
of Am em and B would be somewhat different. However, in that limit Am em 
and B are are negligible compared to m u and anyway, so the shape of the 
curves is little affected.) This gives the condition (1— r)m^ > m u +11.7 MeV. 
If this is satisfied, all neutrons beta decay, even inside nuclei. Since stable 
nuclei having A > 1 cannot be made up only of protons PQ , the only stable 
nucleus in such domains and therefore the only chemical element would be 
1 H. It is unlikely that chemical life can exist in an such "all hydrogen" 
domains. Thus viability requires that one has (1 — r)rrid < m u + 11.7 MeV, 
which rules out the region above the upper green curve in Fig. 4. Another 
constraint comes from the requirement that protons should be stable in 
nuclei. If m u — rrid — m e + Am em > B, then protons even in nuclei would 
be able to decay through the process (A, Z) — ► {A — 1, Z — 1) + n + e + + v e . 
Viability thus requires m u < (1 + r) 777^ + 8. 3 MeV, which rules out the region 
below the lower green curve in Fig. 4. 

A stronger bound comes from the requirement that isolated protons are 
stable against p — ► n + e + + v e . If they were not, there would be no 1 H in 
such domains. Since almost all organic molecules (with a few exceptions, 
such as CO2) contain hydrogen, a no-hydrogen domain would presumably 
have dim prospects for life based on chemistry. This rules out the region 
below the curve 777,^(1 + r) = m u + Am em , shown as the orange curve in Fig. 
4. A slightly stronger bound comes from the condition that the proton in the 
nucleus of hydrogen not convert to a neutron by electron capture p + e~ — > 
n + v e , which rules out the region below the curve m<i(l — r) = m u + Am em , 
shown as the lower red, concave-up curve in Fig. 4. 

If the reaction p+p — * D+e + +v e in the sun is to be exothermic one must 
have 771(7,(1 + r) < m u + Bp + Am em , where Bq = 2.2 MeV is the binding 
energy of the deuteron. A very similar condition comes from requiring that 
the deuteron be stable against the weak decay D — > p + p + e~ + ~P e > namely 
777^(1 — r) < m u + By) + A?7i em . In [1] the binding energy of the deuteron 
was parameterized by 
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B D = 2.2 -a 



( 



m u0 + m d0 



m u + m d 



- 1 



(3) 



where energies are in MeV and the parameter a is very sensitive to the 
model of the internucleon potential. Using two different models, [IJ obtained 
a = 5.5 MeV and a = 1.3 MeV. Using Eq. (3), the condition for the deuteron 
to be stable against weak decay becomes (0.93 — a/ll)irid — (1 — a/ll)m u < 
(3.9 + a). For a = 5.5 MeV, this rules out the region above the upper, 
red, concave-up curve in Fig. 4. As noted in p], it may be possible for 
the hydrogen burning to occur in stars and for nucleosynthesis to proceed 
even if the deuteron is unstable to weak decay, if its lifetime is sufficiently 
long. Thus, a more clear-cut constraint comes from the condition that the 
deuteron be stable against the (much faster) strong decay D — > p + n. The 
condition for this is simply that Bd be positive. Using a = 5.5 MeV, this 
gives m u -\-md < 15.4 MeV, which rules out the region above the upper, red, 
concave-down curve in Fig. 4. 

If the quark masses are too small, then the pion mass becomes small, and 
the nuclear force becomes sufficiently long-range to bind the "diproton" , i.e. 
the 2 He nucleus, as pointed out by Dyson long ago [9]. This would allow 
hydrogen to burn in stars by the reaction p+p — > 2 He, which is very fast as it 
does not involve the weak interactions. Hydrogen-burning stars would have 
lifetimes orders of magnitude shorter than the billions of years presumably 
required for complex organisms to evolve. In [TJ [8] it was estimated that if 
m u + rrid < 0.75(m n o + TUdo) = 8.25 MeV, the diproton is bound. This rules 
out the region below the lower, red, concave-down curve in Fig. 4. 

Altogether, then, we see that the viable region is squeezed down to the 
remarkably small shaded green area bounded by the four red curves in Fig. 4. 
The actual value of (m u ,m<i) observed in our domain is right in the middle 
of it, as shown in Fig. 4. (For perspective, recall that this shaded green 
region is smaller than the dot in Fig. 3, which shows the whole parameter 
space.) 

2.2 Exotic domains with light c quarks 

Exotic possibilities open up if either m c and m s is small. First we shall 
consider the situation where v u is so small that m c becomes comparable to 
or smaller than m,d- (Of course, since the Yukawa coupling ratios are fixed, 
m u becomes about two orders of magnitude less than m^.) We shall call 
these "light-c-quark domains". In such domains the light quarks are d, u 
and c, and there is a baryon octet made up of valence quarks of these three 
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flavors. This octet, shown in Table I, contains baryons of charge 0, +1, and 
+2, but none with negative charges. Therefore these domains are similar to 
ours in the following respects: (1) All nuclei have Z > 0. (2) If a nucleus has 
large A, then its charge Z is also large, since Fermi energy tends to equalize 
the number of baryons of different species. (3) For large Z nuclei, Coulomb 
energy becomes important and gives an upper bound to the size of stable 
nuclei. 



"cc "^c 

(ccd) (cdd) 

3 c + c + £ C +,A+ 

(ecu) (cud) (udd) 

S c ++ P 

(cuu) (uud) 

Table I 



In the light-c-quark part of parameter space one can distinguish three 
regions, which are shown in Fig. 5: (A) For large enough m c , the lightest 
baryons in the octet are those without valence c quarks, namely p and n. For 
nuclei composed of p and n the analysis of the previous subsection applies. 
Consequently, case (A) does not give any anthropically allowed regions of 
parameter space except those we have already discussed. (B) For small 
enough m c , the lightest baryons in the octet are those with only valence 
c and u quarks, namely £+ + and (i.e. uuc and ucc, respectively). 

As we shall see, domains of this sort seem unlikely to be viable. (C) For 
intermediate values of m c , the lightest baryons in the octet are p and 
Some domains of this sort may be viable. 

The sufficient condition for case (A) is that m n < which gives 

< m c + m u + Am' em , where the last term is the electromagnetic contri- 
bution to the splitting of X+ + and n. We may neglect m u in this expression, 
since it is less than a percent of m c . We take Am' em = 4Am em = 6.8 MeV. 
(Note that the observed splitting between S+ + and £° is 0.3 MeV, so that 
if the quark mass contribution is 2(m u — ma) = —6 MeV, it implies an 
electromagnetic contribution of 6.3 MeV. Of course, the wavefunctions for 
the charmed baryons would be different in the light-c-quark domains.) The 
condition is then m c > 2m^ — 6.8 MeV. This is the region to the right of 
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x d =ln (m d /M P1 ) 




x u = ln (m u /Mpi) 



Figure 5: The "light-c-quark" region of parameter space. The green shaded 
region may be viable. The dot is our domain. 



the lower red curve in Fig. 5. In this region, the analysis we gave in section 
2 applies, and one expects that only anthropically allowed region is that 
shown in Fig. 4, which includes our domain (indicated by the blue dot in 
Fig. 5). 

The condition for case (B) is that m~++ < m v , which gives 2m c + 
3Am em < m<i + m u , or approximately m c < — 2.5 MeV. This is 

the region to the left of the upper red curve in Fig. 5. This region seems un- 
likely to be viable since all stable nuclei would have Z > 2. (It is easily seen 
that an isolated p will beta decay to as will a p bound with neutral 

baryons.) Consequently, there would be no chemical analogue of hydrogen. 

The most interesting case is (C), which lies in between the two red curves 
in Fig. 5. Here the lightest two baryons are p and In such a domain, 



13 



the only Z = 1 nucleus consists of a single proton, i.e. H. In order for 
hydrogen to be stable, it is required that the decay p — > £+ + + e~ +V e be 
energetically forbidden, which gives the condition m c > m^l — r) — 3Am em . 
This rules out the region above the blue, concave-up curve in Fig. 5. A 
second condition is that diprotons not be bound, as otherwise hydrogen- 
burning stars would burn too fast, as noted above. This gives the bound 
that rrid + m u > 8.25 MeV, as given previously. This rules out the region 
below the nearly horizontal blue line that barely misses our domain in Fig. 
5. 

A third condition is that the A = 2 nucleus composed of a p and a 
S+ + be bound. Let us call this analogue of the deuteron D. The potential 
between the p and the S+ + at large distances is given by one-D-meson 
exchange rather than one-pion-exchange, and is thus controlled by m c + 
rather than by m u + m^. Using Eq. (3), the binding energy of D can be 
written = 2.2 — E c — a[(m c + m,d)/(m U Q + rrido) ~ l]j where E c is the 
Coulomb interaction energy of the baryons in the D. The condition that 
D be bound is then m c + < 11 [(2.2 — E c )/a + 1), where energies are in 
MeV. If a = 5.5 MeV and E c = 2 MeV, this rules out the region above the 
blue, concave-down curve in Fig. 5. 

The combination of the three constraints leaves only the very small tri- 
angular region shaded green in Fig. 5. In fact, since the three curves depend 
on quantities that are not well known, such as a, E c , m U Q, m^ , it is not 
really clear that any area is actually enclosed by them. Furthermore, even 
if there is, it does not necessarily mean that this area is an anthropically 
allowed region of parameter space. It only means that three of the more ob- 
vious disasters for viability do not happen. One other potential problem for 
viability is that nucleosynthesis may be inhibited in such domains because of 
large Coulomb barriers, due to the fact that nuclei are composed of charge- 1 
and charge-2 baryons. A much more serious potential problem for viability 
is that Gf is about four orders of magnitude larger in this region of pa- 
rameter space than it is in our domain. This would make the ll pp reaction" 
p + p^D + e~+V e proceed in hydrogen-burning stars very much faster 
than the ordinary pp reaction does in our domain. That would presumably 
be disastrous for the same reason that the diproton being bound would be. 

As we shall see in section 3, where we discuss probability distributions 
in parameter space under some simple and reasonable assumptions, these 
potentially viable light-c-quark domains occur less commonly in the universe 
than domains like our own. Therefore, unless a typical domain of this sort 
is more viable than ours, the existence of such domains should not affect 
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the anthropic explanation of the observed values of m u and rrid that we are 
exploring here. 

2.3 Exotic domains with light s quarks 

We now consider the situation where Vd is so small that the mass of the s 
quark is comparable to or smaller than m u . (Of course, then rrid <C m u , since 
we are assuming that the ratios of Yukawa coupling are fixed.) We will call 
these "light-s-quark domains". In such domains, the lightest three quarks 
are d, s, and u, as in our domain, but the corresponding flavor SU (3) is a 
very good symmetry, since m s is small. Consequently, there is the familiar 
baryon octet, shown in Table II, but it is more nearly degenerate than in 
our domain. 

H° £+ 

(uss) (uus) 

a- £°,A° p 

(dss) (uds) (uud) 

£- n 

(dds) (udd) 

Table II 



There are three cases to consider for light-s-quark domains (just as there 
were for light-c-quark domains). These are shown in Fig. 6. (A) For large 
enough m s , the lightest two baryons are the two non-strange baryons p and 
n, as in our domain. The condition for this is that m s - > m p , which gives 
m s + 171(1 > 2m M . This is the region above the upper red curve in Fig. 6. (B) 
For small enough m s , u is heavier than both s and d and the lightest two 
baryons are £~ and H - , which have no valence u quarks. The condition for 
case (B) is that m s - < m n , which gives m u + > 2m s + Am em . This is 
the region below the lower red curve in Fig. 6. (C) For intermediate m s , 
the lightest two baryons would be n and E - . This is the region between the 
two red curves in Fig. 6. 

For case (A), where the lightest baryons are p and n, as in our do- 
main, the analysis of viability we presented earlier applies, and so the only 
anthropically allowed region is that shown in Fig. 4, which includes our do- 
main (indicated by the blue dot in Fig. 6). Case (C) is interesting because 
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Figure 6: The "light-s-quark region of parameter space. Some of the region 
below the green line and to the right of the rightmost red curve may be 
viable. The dot is our domain. 



it might lead to domains in which the S~ plays the role that the proton 
plays in our domain. For example, a "hydrogen atom" in such a domain 
would consist of a T,~ and an e + . (One would expect antileptons to pre- 
dominate in such domains because of charge neutrality.) There are at least 
three constraints that m u and m s must satisfy in case (C). 

First, the S~ in the nucleus of the "hydrogen" atom must be stable 
against beta decay into n + e~ + V e . This gives the constraint m u > (1 — 
0.05r)m s + Am em . This rules out the region to the left of the blue curve 
that is concave to the right in Fig. 6. 
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The second constraint is that the n and S _ must be able to bind to 
make the analogue of the deuteron. The long-distance part of the potential 
between these nucleons is controlled by the mass of the kaon, and thus by 
m u + m s . Thus the condition for the n and S~ to bind is m u + m s < 15.4 
MeV, if a = 5.5 MeV. This rules out the region above the blue curve that 
is concave to the left in Fig. 6. The third constraint is that the di-X~ not 
be bound (the analogue of the diproton not being bound in our domain). 
This gives m s + m,j > 8.25 MeV, which rules out the region above the nearly 
horizontal blue line in Fig. 6. Altogether, there is no allowed region left for 
case (C). However, this depends on the value of a, and other parameters that 
are not precisely known. If a is less than about 3 MeV, there may be a small 
region that satisfies the three constraints. Of course, even if there were, that 
hardly proves that such a region is anthropically allowed; it merely means 
that some obvious disasters for viability are avoided. 

We turn now to case (B). In this case the lightest baryons are T,~ and 
S _ , and E~ is the lighter of the two. The only obvious constraint is that 
the analogue of deuterium (i.e. a bound state of S~ and E~) be allowed to 
form. For a = 5.5 MeV, this gives the constraint m s + = 1.05m s < 15.4 
MeV. This rules out the region above the horizontal green curve in Fig. 6. 

There may be some anthropic upper bound on m u , coming from the fact 
that very large m u corresponds to very small Gf; however, that is not clear. 
If v u is superlarge and therefore Gf ~ v~ 2 is supersmall, rather strange 
situations are possible. Even though Gf is supersmall, the decay rates of 
u, c, and t quarks, being of order G 2 F m\ c t ~ v u will be very fast, and the 
only quarks around will be d, s, and b. The amplitudes for transitions of d, 
s, and b into each other will be suppressed by a factor of v~ 2 and therefore, 
in effect, there will be separate conservation of (i-number, s-number, and 
6-number. (Sphaleron processes will, of course also be highly suppressed.) 
The story is similar for the leptons; there will be separate conservation of 
L ei L^, and L T . 

In this case, much depends on the details of primordial lepto/baryogenesis. 
It is possible that lepto/baryogenesis leads to asymmetries for d, s, and b 
quarks that differ in sign. For example, there could be an excess of d over 
d, but of s over s and b over b. In that case, "ordinary" hadronic matter 
would be made up of d, s, and b. In particular, there would be stable K°, 
B° mesons, £~ baryons, and positively charged antibaryons consisting of 
valence s and b. These could capture leptons of positive or negative charge 
(as the case may be) and make atoms. In the same way, there could be 
asymmetries of different sign for different flavors of lepton. For instance, 
there could be an excess of e~ and (i~ but of r + . Thus there could be atoms 
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consisting of an e~ orbiting a r + . Clearly, all sorts of rich and strange pos- 
sibilities for chemistry might exist. And if the chemistry is rich, then it is 
plausible that life might be possible. However, in order to have any nuclei 
with Z > 1, one still needs two baryons to be able to bind. There is still, 
therefore, an anthropic upper bound on + m s for case (B), as discussed 
above. As is clear from the shape of the probability contours shown in Fig. 
3 and derived in section 3, this will put an upper bound on the probability 
density of viable small-s-quark domains. 

In sum, there are regions of parameter space with m s of order a few MeV, 
and in particular case (B), where we have found no very obvious convincing 
argument against viability. There may also be a small anthropically allowed 
region where m c is of order a few MeV. As we will see in section 3, for certain 
plausible assumptions about the probability distribution in the lnm u — In 
plane, the possibility of such exotic viable light-s-quark or light-c-quark 
domains would not affect the anthropic explanation of the observed values 
of m u and being explored in this paper. 

2.4 The case of> 2 > 

In domains where fj, 2 > the vacuum expectation value of H < is of order 
/^//u 2 and therefore extremely small compared to its value in our domain. 
This means, of course, that in these domains the quark and lepton masses 
are correspondingly small. (For example, if fj, 2 = +(100GeV) 2 and tan/3 is 
the same as in our domain, then the electron mass is only about 5 x 10~ 4 eV.) 
This has important consequences for life based on chemistry. The reason is 
that the energies of typical atomic and molecular transitions are of order, 
or smaller than, a 2 m e , and therefore so are the temperatures needed for 
biological processes to occur. (That is why in our domain the temperature 
needed for life is about 300 K = 10~ 3 a 2 m e o.) It follows that before life 
based on chemistry can appear the universe has to have cooled down to a 
temperature much smaller than a 2 m e . 

Because m e is so small in positive-^ 2 domains, it was argued in [1] that 
by the time the universe cools sufficiently for chemistry-based life to exist 
in those domains various disasters might already have occurred, such as all 
baryons having decayed away or all stars having burned out. Those ar- 
guments depended crucially on the assumption that universe expands as a 
power of t. However, it now appears that after dark energy comes to dom- 
inate the universe will expand exponentially in t. Once the expansion be- 
comes exponential, it does not take long for the universe to reach extremely 
small temperatures. For example, consider a domain where m e = 5 x 10~ 4 
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eV. If there were no dark energy, it would take such a domain about 10 
years to reach temperatures of order 10 _3 a 2 m e . On the other hand, if there 
were the observed amount of dark energy, it would take only about 10 11 
years. The difference is even more dramatic if p 2 is larger. In a typical 
domain where p? = (10 10 GeV) 2 , it would take about 10 44 years to reach 
biochemical temperatures if there were no dark energy, but only about 10 13 
years if there is the amount of dark energy we observe. 

The existence of dark energy means that the specific astrophysical and 
cosmological anthropic arguments used for the positive-/^ 2 case in pQ are not 
valid. However, there are other anthropic arguments based on astrophysical 
and cosmological considerations that become possible precisely because of 
the existence of dark energy, as we will now see. For the rest of the paper 
we will assume that the dark energy has w = — 1 and will refer to it as "the 
cosmological constant" . Our arguments should not be qualitatively affected 
if w is not exactly —1. 

It seems reasonable to suppose that life requires the existence of planets 
upon which to evolve. A planet is supported against gravitational collapse 
by the electrons, and has therefore a typical density that goes as p p i a net ~ 
Am p (am e /2) 5 , where A is the average number of baryons per nucleus. A 
planet will be ripped apart by the inflation caused the cosmological constant, 
unless 

PA < Ppianet ~ Aa 3 m 3 e m p . (4) 

If pa is truly a constant of nature, so that it the same in every domain, then 
its measured value of about 10 -123 Mp f means that planets can exist only if 
m e > 6 x 10 -5 eV. For tan (3 having the same value as in our domain, this 
means that p 2 has to be less than about (300 GeV) 2 . This is significantly 
stronger than the bound given in Eq. (6) of [I] (a bound that came from 
not having baryons decay away before the universe cooled to temperatures 
where chemical life could exist). 

However, it is more interesting to suppose that p\, as well as the Higgs 
mass parameters, varies among domains. Indeed, this comports with the 
idea that all the dimensionful parameters of the Lagrangian are determined 
anthropically. In that case one has to consider simultaneously the proba- 
bility distributions of p 2 , tan/3 and p\. This will be done in section 4. It 
will be seen that the viability of a significant part of the positive p? pa- 
rameter space is left unaffected by the "planet argument". However, most 
of that part of parameter space can be anthropically ruled out by a second 
argument, which is based on structure formation. 
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That second argument assumes that structure cannot begin to form until 
recombination occurs. If by the time of recombination the density of the 
universe is already dominated by the cosmological constant, then structure 
will never form. The recombination temperature T rec is about 0.05a 2 m e , so 
the constraint is that structure cannot form unless 

PA < PB{Uec) ~ m p r] B T^ ec ~ W~ 14 a e m 3 e m p . (5) 

This appears to be a much more stringent constraint than the one in Eq. 

(4) . However, it is not clear that it can be used for m e smaller than about 20 
cV. The point is that, if m e < 20 eV, then when matter begins to dominate 
over radiation (at T ~ 1 eV) the electrons will still be almost relativistic, 
and the thermal number density of electrons and positrons will be larger 
than that of the baryons by a factor of about 10 10 e _me / T . There will then 
be of order 10 10 electrons and positrons for each baryon. The charges of the 
baryons will be Debye screened, and it would seem that baryons can already 
begin to condense at that point, since Sps/pB would not have a significant 
effect on the density of the leptons and photons. We will therefore only use 
Eq. (5) when m e > 20 eV. The anthropic bounds coming from Eqs. (4) and 

(5) will be analyzed in section 4. 




3 Probability distributions in the In m u — In rrid plane 

Since anthropic explanations of relations among parameters are based on 
probabilities, one must be able to make reasonable hypotheses about prob- 
ability distributions within the parameter space or "landscape". The pa- 
rameters that we are assuming to vary among domains of the universe are 
those in the Higgs mass-squared matrix Mfj. This matrix acts in the space 

It is convenient to parameterize this matrix as follows 

M 2 = ( ao + a3 ai ~ m2 ). (6) 
y ai + ia 2 cio- a 3 J 

First, we will assume that the four real parameters a a (which have di- 
mensions of mass squared) are "naturally" superlarge (either 0(Mq UT ) or 
0(Mp £ )) and of the same order. Call their natural scale M 2 . The simplest 
assumption is that the four a a are independent and that each has a flat 
probability distribution in the interval [— M 2 ,+M 2 \. (However, symmetry 
may somewhat suppress the off-diagonal elements of Mfy We will discuss 
the effects of this later.) We will take oq to be positive. (If it is negative and 
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very large the weak interactions are broken at very large scales, whereas if 
it is negative and small then breaking the weak interactions at small scales 
would require that 03 also be small, which is an extra fine tuning.) It is 
convenient to think of the parameters a±, 02, and 03 as forming the compo- 
nents of a vector a in a three-dimensional space. As these three parameters 
have been assumed all to have flat probability distributions in the inter- 
val [— M 2 , +M 2 ], the probability distribution is rotationally invariant in the 
space of a. 

We shall denote the larger eigenvalue of M 2 - by M 2 and the smaller 
eigenvalue by fi 2 . (By our assumption about ao, M 2 is positive.) Since M 2 
is generally much larger than /j, 2 , the heavier Higgs doublet decouples; so the 
parameter /j 2 is in effect just the /j, 2 of the Standard Model and controls the 
breaking of the weak interactions. Because the cases of anthropic interest 
are those where \/j, 2 \ <C M 2 , we may write 

2 ^ detM 2 _ a 2 -a 2 
^ ~ tvM 2 ~ 2a ' 1 j 



21 M 2 , 1 \/«?+«9 

tan 2/3 = ' = ^ = tan 9 a , (8) 

M 2 U - M 2 d a 3 

where (3 is the angle between the mass basis of the Higgs doublets (H > ,H < ) 
and the original flavor basis (H u , H^) (i.e. H < = e %a sin j3H u + cos (3Hd and 
tan/3 = ^ ), and 9 a is the angle between a and the 03 axis. There are two 
cases to consider, fi 2 < and ^u 2 > 0. 

3.1 Probability distribution in the case of fi 2 < 

In the case of negative /j, 2 , the vacuum expectation value that breaks the 
weak interactions is v = ^\fi 2 \/X = (2A) -1 / 2 |(oq — a 2 )/ao| 1 ^ 2 , where a 
is the magnitude of a, and A is the quartic self-coupling of the lighter 
Higgs doublet. Then we may write m u = Y u v u = Y u v sin/3 = (Y U /V2X) 

sin If^J (a 2 , — a 2 )/ao, and similarly nid = Y^va = Yjv cos (3 = (Yd/ V2X) cos 

\J (oq — a 2 )/ao- Defining K = Y U /V2X and recalling the definition of k in 
Eq. (2), one has 



m u = Ksm^f^J(al - a 2 )/a , 
m d /k = Md~ = K cos?fJ(al - a 2 )/a . 
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Since, by our assumptions, the probability distribution is rotationally invari- 
ant in the space of the vector a, one may write the probability distribution 
in spherical coordinates (a,9 a ,4> a ) in that space: 



P(ao, a, & a , 4>a) dao da d6 a d(f> a = N dao a 2 da sin 6 a d6 a d(f> a , 



(10) 



where TV" is the appropriate normalization factor. We would like to compute 
the probability distribution in the space of m u and m d . Thus we write 
P(ao,a,9 a )dao da d6 a d(p a = P(ao,m u ,rrT d )daQ dm u drfQ d(p a , and compute 
the Jacobian for the transformation of coordinates: 



J 



d(ao,m u ,m d ) 



d(a ,a,9 a ) 



l -K\a/a ). 



This yields the result 

P(ao,m u ,rnd) = 2N K~ 2 ao a sin# a . 



^From Eq. (9), one has a = y Oq — a§K~ 2 (m\ + m/) and 

2m u m2 



sin On = 



m\ + m d 2 ' 



This yields 



(11) 



(12) 



(13) 



P(a ,m u ,m d ) = AN K a v/ag - a Q K 2 {m 2 + m d 2 ) u , . (14) 

v + m d z 

For the regions of parameter space that are anthropically interesting (i.e. 
where life is likely to be possible) one has K -2 m^, K^m^ 2 <C ao ~ M 2 , so 
that the previous equation can be written P(ao, m u , rn^) = \N K~ 2 a^m^fn^j {m 2 l + 
m^j 2 ). One then simply integrates over the variable ao to get 

T3 ( x , , ,m u m- d dm u dmr d 2 

P{fn u , m d ) dm u drag = N ^ == ^ — , for /i < 0. (15) 

mi + m d z 

It is more interesting to express the probabilities in terms of the loga- 
rithms of the quark masses. Defining x u = ln(m u /M*), x d = ln(m d /M*), 
and a?d = ln(TOd/M*) = x d — ln/c, one has 

P(x u , x d ) dx u dx d = N" -^r &m dx u dx^, for fi 2 < 0. (16) 
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The contours of constant probability given by this expression are plotted 
as the concave- up curves in Fig. 3. Note the shift represented by the fact 
that P depends on Td = Xd — Ink, where k ~ 10 2 . This shift means that 
domains where = m^o and k~ 1 m U Q < m-a < ^uo are just as likely to 
occur as domains like ours. The reason for this is simple. In our domain, 
one has rrid ~ m u , even though Yd = kY u » Y u . This means that v u /vd = 

tan/3 ~ k S> 1, so that v = ^Jv 2 + v \ = v u . Reducing m u from its value in 
our domain (m u o) while holding rrid fixed, means reducing v u while holding 
Vd fixed. This has the effect of reducing v and thus /j, 2 , which increases the 
fine tuning. However, it also has the effect of making the ratio v u /vd closer 
to one, which reduces the fine-tuning. (The factor of s'm9 a = sin 2(3 in P 
prefers equal values of v u and Vd-) These effects cancel, so that it does not 
"cost probability" to move to the left in Fig. 2 from the point representing 
our domain down to values m u ~ k~ 1 m U Q. Decreasing m u more than that 
does cost probability, however, because it means making v and fi 2 more fine- 
tuned without any compensating effect from tan/3, which indeed starts to 
move away again from 1 towards 0. 

Thus, to explain why we don't observe a value of m u as small as k~ 1 m U Q ~ 
10~ 2 m u o, one cannot argue that such domains are less common. Rather, 
one has to argue that they are less viable. In fact, we have presented such 
arguments in section 2. What makes the region of parameter space with 
m u <C m U Q less viable is a combination of constraints coming from reactions 
in hydrogen-burning stars like the sun, as one sees from Fig. 4. If m u is even 
slightly smaller than m U Q, it either makes the pp reaction p+p — > D + e + + v e 
endothermic (since protons become lighter relative to neutrons), or it makes 
the diproton a bound state, allowing the reaction p + p — > 2 He, or it has 
both effects. Another effect that is probably deleterious to the chances of 
life is that as v u gets smaller v gets smaller, as mentioned, and so Gf gets 
larger. This would make the pp reaction go much faster (since it involves the 
weak interaction), presumably making stars like the sun burn much faster. 
However, we have not calculated this effect, since the reaction p + p — > 2 He 
being allowed is an even more drastic effect. 

3.2 Probability distribution in the case of ji 2 > 

If fi 2 is positive, then the light Higgs doublet H < only acquires a vacuum 
expectation value due to quark-antiquark condensation at /„-, the chiral- 
symmetry-breaking scale of QCD. Specifically, the lighter Higgs doublet has 
the coupling H < (Yt sin (3tt + cos (3bb). After QCD chiral symmetry break- 
ing this produces a linear term in the effective potential for H < that is of 
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the form (Y t sin/3 + Y& cos f3)f 3 H < . Since we have assumed that Yt ~ Yf, ~ 1, 
the coefficient {Yt sm P + ^5 cos /3) is a slowly varying function of /3 that is 
of order 1. We will denote this j(/3). 
Then instead of Eq. (9) one has 



m u = Y ul {f3)fl sin % ((a 2 - a 2 )/2a ) 



(17) 



m d /A;= W^ = y u7 (/?)/3 cos ^(( a 2_ a 2) /2ao ) 
Thus, the Jacobian is given not by Eq. (11) but by 

J = (Y ul {0)f n )Hava 2 /{al - a 2 ) 3 . (18) 

There may be values of f3 for which 7(/3) <C 1 due to an almost exact 
cancelation between the Yt term and the Yj, term. However, that simply 
makes m u and m d small. Since smaller values of m u and m d are the the 
most probable ones anyway in the case of /x 2 > 0, the possibility of such 
cancelations has virtually no effect on the probability distribution we shall 
derive. Thus we will get qualitatively the right behavior for P(m u , m^) if we 
treat 7 as a constant, which we will now do for simplicity. Then following 
the same steps as before one ends up with 

~ .m u m^dm u drrTd 2 

P{m u , m d ) dm u drrQ = N - — 2 , , for fi > 0. (19) 

(7774 + md ) 1 

and, in terms of the variables x u and x^, 



~ 6 -(-2 

P(x u ,x d ) dx u dx d = N" ^ - e2 ^ 5/2 dx u dxZ, for fj 2 > 0. (20) 

The probability contours for both positive and negative /j? that we have 
just derived have been displayed in Fig. 3. For negative /i 2 the probabilities 
increase upward and to the left (i.e. towards larger m u and m d ), whereas for 
positive /i 2 the probabilities increase downward and to the left (i.e. toward 
smaller m u and m d )- The directions of increasing probability are shown 
by the arrows and the contours with equal probability to our own domain 
are shown by the red curves. We see from this that the potentially viable 
"light-s quark domains" and "light-c-quark domains" are less common than 
domains like our own. 

In deriving the probability distributions given in Eqs. (16) and (20), we 
assumed that all the elements of the Higgs mass matrix M 4 2 have the same 
natural scale M 2 . However, the "natural flavor conservation" [5] pattern of 
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Yukawa couplings that we have assumed in Eq. (1) (whereby H u couples to 
up-type quarks and couples to down-type quarks and charged leptons) 
is not "technically natural" unless there is a symmetry K under which H u 
and Hd transform differently. The off-diagonal elements of M 2 j necessarily 
break this symmetry. If K is a global symmetry, then one can simply regard 
the d = 2 term M 2 d HlH d as a soft breaking of K. In that case, the pattern 
of Yukawa couplings in Eq. (1) is still technically natural. It would be 
reasonable to assume that the natural scale of M 2 d was the same as that of 
MiandMf,. 

On the other hand, one might imagine that K is spontaneously broken. 
For example, suppose M 2 d arises from a X-invariant term eH^H^S 2 , where 
S has a potential V(S) = \Xs(\S\ 2 ) 2 + fi^\S\ 2 . Then it would be consistent 
with our approach to say that ji 2 s scanned and had a natural scale M 2 . In 
this scenario, there are two equally probable cases, fig > and fi 2 s < 0. For 
^5 < 0, the natural scale of M 2 d is e\g M 2 . If eAg is small, the probability 
distribution P(x u ,Xd) is suppressed when v u and are comparable (more 
precisely, when they are within a factor of eA^ 1 of each other), i.e. near 
the line x u = x~d in Fig. 3. This would not qualitatively affect our analysis 
unless eA^ 1 were very small (less than about 10~ 2 ), in which case the poten- 
tially viable small-c- quark and small-s-quark regions would start to have 
probability densities rivaling that of our domain. If jj? s > 0, then M 2 d = 
and quite a different situation results. A single fine-tuning to make the weak 
scale small would leave either M 2 U or M dd of order +M 2 . If M 2 U ~ +M 2 , 
the u, c, and t quarks would have negligible masses, while the down quark 
masses would be set by the scale of weak interaction breaking. We have not 
investigated the viability of such domains. 

In our discussions of probability distributions, we have not considered 
the renormalization group running of A and fj, 2 from the scale M* down to 
the scale v. This would affect the probability distributions by factors of 
order 1. However, it should not qualitatively affect our discussion. 



4 Simultaneous tuning of the Higgs mass matrix 
and A 

In section 2.4, two anthropic constraints on the relation between the mass of 
the electron and the cosmological constant were discussed. These eliminate 
as non-viable the great majority of the positive-^ 2 parameter space. How- 
ever, in order to express these as anthropic constraints on the Higgs mass 
parameters (or equivalently on x u = ln(m u /M*) and Xd = ]n(md/M*)) one 
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must do an analysis that simultaneously takes into account the probability 
distributions of x u , x d and p\. That is done in this section. Both of the 
anthropic constraints on p\ discussed in section 2.4 and given in Eqs. (4) 
and (5) require p\ to be less than bounds that are proportional to rv? e and 
thus to e 3xd . If we assume that the probability distribution of p^ among 
domains is flat in the interval [—Mp^+Mp^, then the probability distri- 
bution given in Eq. (20) gets modified by a changed normalization factor 
and by an extra factor of e 3xd to account for the tuning of the cosmological 
constant required to satisfy Eq. (4) or (5): 

P'(x u , x d ) dx u dx d oc — — - e2 ^ b/2 e 3xd dx u dxl. (21) 

The prime on this probability density indicates that this is a conditional 
probability: it is the probability density for being at (x u ,Xd) subject to the 
condition that p\ has been tuned to a small enough value to satisfy either 
Eq. (4) or Eq. (5). If the condition is that Eq. (4) be satisfied, we will 
call this Ppi anet , if the condition is that Eq. (5) is satisfied, we will call it 
P' 

structure- 

The constant of proportionality in Eq. (21) depends on which condi- 
tion is being enforced. To determine the constants of proportionality, it is 
convenient to start by finding the points in the positive-^ 2 and negative- 
p 2 regions of the x u — Xd plane where x u = Xd and where the Higgs mass 
matrix M 2 - is fine-tuned to the same degree as in our domain (without any 
condition on p\). We will call these point S+ and S-. (The S stands for 
symmetric, since x u = x d there, and the subscript tells whether it lies in the 
p? > or p? < region.) These points are shown in Fig. 7, where the point 
representing our domain is labeled D. Suppose that Eq. (5) implies that at 
S+ the cosmological constant has to satisfy p\(S + ) < V structure PAO, where 
PAO is the value observed in our domain, i.e. that p\ must be fine-tuned at 
S+ to be less than it is in our domain by a factor smaller than "p structure- 
On the other hand, the Higgs mass parameters are equally fine-tuned at 
5+ and in our domain, by definition of S±. So the conditional probability 
density P' s tructure a t S+ divided by the conditional probability density P' in 
our domain (where p\ is tuned to pao) is just given by p s tructure- Therefore, 
if we define the shifted variables x u = x u — x u (S + ) and x~d = Xd — x d (S + ), 
we may write 

pl (V o5/2 2x u +5x d 

r structured u, -^d) * c / 99 n 

p/( D) Pstrueturc ^ + ^5/2 ' ^> 
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Figure 7: The i? > region (approximately v u ,Vd < f n ) is below and to 
the left of the blue curve. The arguments of section 2.4 rule out all but the 
shaded blue region. 



Note that at S+, where x u = x~d = 0, the right-hand side of the equation is 
just p structure- In order to compute the coefficient p s tructure, it is necessary 
to determine the value of m e (S+) and substitute it into Eq. (5). 

At the point D, i.e. in our domain, one has v u = v u q ~ 174 GeV, 
Vd = Vdo — 4 GeV (if we assume Yy, ~ Y t ~ 1), and m e = m e Q = 
Y e v d0 ~ Y e (4 GeV). Since P(S+) = P(S-), Eq. (16) implies that v u (S-) = 
Vd(S-) ~ V2v d0 ^ 4\/2 GeV, and therefore v(S-) ~ 8 GeV; so that 
m e (5_) = (Y e /V2)v(S-) ~ Y e (iV2 GeV) ~ V2rn e0 = 0.7 MeV. On the 
other hand, at the point S + , the electron mass is given by m e (5 + ) = 
Q eA/2) {fl /n 2 { S+)). But by definition, /j 2 {S + ) = /u 2 (5_) and therefore 
^fi' 2 (S+)/X = y/tf(S-)/\ = v(S-) ~ 8 GeV. Taking A ~ 1 we have 
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m e (S + ) ~ m e (S'_)(/ 7r /(8GeV)) 3 ~ 1 eV. Substituting this into Eq. (5), 
one has p\(S + ) < 10 _7 pao- In other words, 

Pstructure = 10 . (23) 

The constraint on p\ coming from Eq. (5) only applies if m e > 20 eV. 
The region with m e > 20 eV lies above the horizontal straight red line in 
Fig. 7. The region of positive p? corresponds to the area below and to the 
left of the blue curve in Fig. 7. (This curve arises from the fact that for very 
small positive p 2 , the value of (H<) is set by the competition between the 
quartic term (A/4)i7< and the linear term for H < coming from the QCD 



condensates of qq. This gives a maximum value for y v u + v d = (i?<) — f n 
in the positive-// 2 region.) It is easy to see from Eqs. (22) and (23) that the 
entire region satisfying both /i 2 > and m e > 20 eV (shaded gray in Fig. 
7) is much more fine-tuned than our domain, i.e. Pstructure/ P '(D) *C 1. 

Turning now to the condition that planets be able to form (Eq. (4)), the 
calculation is completely parallel to the one just done. To find p p i an et on e 
substitutes m e (S + ) into Eq. (4) and divides by /?ao> which gives 

Ppt anet (Xu,X d ) 

P'{D) <Pplanet ( e 2x u +e 2x d ^5/2- ^> 

and 



Pplanet ^ 3 X 10" A, (25) 

where A is the average atomic weight of a nucleus in the planet. It is easy to 
see that the region to the left of the vertical orange line or below the slanted 
orange line are more fine-tuned than our domain, i.e. P' p i ane t/P'{D) < 1. 
The shape of the orange line arises from the fact that when £d > x~ u , the 
expression in Eq. (24) is approximately proportional to simply e 2x " , whereas 
when x u > x^, the expression in Eq. (24) is approximately proportional to 

The two arguments given in section 2.4 have excluded (as being more 
fine-tuned than our domain) the whole of the positive-// 2 region except for 
the region shaded light blue in Fig. 7. However, that is not good enough. 
It is necessary to exclude the entire positive-// 2 region either as not viable 
or as more improbable than our domain, if the anthropic explanation of the 
observed values of m u and rrid is to be tenable. Other arguments that were 
suggested in [1] hold some promise of being able to do this. For example, 
as noted there, it might be that the peculiarities of nuclear physics in the 
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positive-^ 2 region may lead to runaway nucleosynthesis. However, the situ- 
ation for fi 2 > has many unclear aspects, and more thought must be given 
to it. 

5 Conclusions 

We have investigated the possibility that the entire mass matrix of the Higgs 
fields in a two-Higgs-doublet model varies among domains of a "multiverse" . 
This means that both the weak scale v = ^Jv 2 + v\ and tan j3 = v u /vd vary 
among domains. Making plausible assumptions about probability distribu- 
tions of parameters in the "landscape" and about the requirements for life, 
we have found that almost all of the In m u — In space is either anthrop- 
ically non-viable or much more fine-tuned (and therefore more rare in the 
multiverse) than our domain. The potentially viable regions are islands in 
parameter space. One of these islands is ours, and it is a very small island 
indeed, as seen in Figs. 3 and 4. If this is the only island of viability, there 
is the possibility of accounting for two otherwise not easily understood facts 
about the world we observe, namely that the weak scale is comparable to 
the strong scale, and that m u /m^ ~ 1. 

Whether this is in fact the only island of viability is not yet clear. There 
are possible islands of viability in the "light-c-quark" and "light-s-quark" 
regions. However, even if islands of viability do exist in those parts of 
parameter space they are less probable in the landscape than the island of 
viability in which our domain is located, according to the analysis of section 
3. More significant is the possibility of an island of viability in the /u 2 > 
region of parameter space. 

An interesting new twist has been given to the discussion of the positive- 
/j? case by the discovery of dark matter /cosmological constant. This discov- 
ery ties together the anthropic analysis of the Higgs mass given in [I] and 
[8] and the anthropic analysis of the cosmological constant given in [TO] . We 
have seen that some of the /u 2 > constraints are invalidated by the existence 
of dark energy, but other and equally strong constraints that depend on the 
existence of dark energy come into play. However, the constraints we have 
discussed are not by themselves enough to rule out the entire /i 2 > region 
as non-viable. This seems to be the main loophole at present in the attempt 
to explain m u ~ rrid and the value of the weak scale anthropically. However, 
it seems probable that much stronger arguments exist for the case of /i 2 > 0, 
and that the entire positive-// 2 region of parameter space is non-viable. 

We have assumed that only the dimensionful parameters of the La- 
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grangian "scan" , i.e. vary in the multiverse, and that only they are anthrop- 
ically tuned. This seems like a comparatively conservative assumption, since 
the dimensionful parameters of our current theories are the least well under- 
stood, are the most (apparently) fine-tuned, and are the most amenable to 
anthropic explanation. However, obviously, other assumptions are possible. 
The question for us is whether there exists at least some set of plausible as- 
sumptions under which the observed relation m u ~ md can be anthropically 
explained. 
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